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' Abstract 
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Motivated by recent experiments with superradiant Bose-Einstein Condensate (BEC) we con- 
sider simple microscopic models describing rigorously the interference of the two cooperative 
. ^ | phenomena, BEC and radiation, in thermodynamic equilibrium. Our resuts in equilibrium con- 
firm the presence of the observed superradiant light scattering from BEC: (a) the equilibrium 
superradiance exists only below a certain transition temperature; (b) there is superradiance 
and matter-wave (BEC) enhancement due to the coherent coupling between light and matter. 



Keywords: Bose-Einstein Condensation, Superradiance 
PACS: 05.30.Jp, 03.75.Fi, 67.40.-w. 

Running title : BEC and superradiance 



1 Email : Joe.pule@ucd.ie 

2 Email : andre.verbeure@fys.kuleuven.ac.be 

3 Email : zagrebnov@cpt.univ-mrs.fr 



1 



1. This letter is motivated by the interest in Bose-Einstein Condensation (BEC) of bosons in 
traps and in particular, by the recent discovery of the Dicke superradiance and BEC matter- 
wave amplification, see e.g. [T], [2], In all these experiments, the condensate was illuminated 
with a Q-mode laser beam (so-called "dressing beam" ) and then the BEC atoms scatter photons 
from this beam into another mode and receive the corresponding recoil momentum producing 
a coherent "four- wave mixing" of light and atoms 3J. Notice that Girardeau in 1978, [TT], had 
already anticipated the possible existence of these phenomena. 

The irradiation of a Bose-Einstein condensate can be considered as an external action on it. 
The effect of other external agents have been studied, for instance the influence of boundary 
conditions (jlj and and that of an external field, scaled with the volume so as to retain the 
space homogeneity [H] . In the first case it is well known that attractive boundary conditions en- 
hance condensation. Here condensation occurs even in one dimension and is a consequence of a 
discrete point in the spectrum. In the other case one also finds an enhancement of condensation 
compared to that of the free Bose gas depending on the behaviour of the potential. 

In the present letter we consider two simple models which show that some particular inter- 
actions with a one-mode radiation can enhance the conventional BEC of the perfect Bose gas 
(PBG). This can be interpreted as a coherent coexistence of the BEC and the condensation of 
photons, which is a type of equilibrium Dicke superradiance induced by BEC, see e.g. j7], [S], 

M- 

We consider a system of non-interacting bosons (PBG) of mass m enclosed in a cubic box 
A C ~R d of volume V — |A| = L d centered at the origin with periodic boundary conditions. The 
Hamiltonian for the PBG is 

T A := ^ e kO*k a k (1) 

where 

{2irn 
keR d : k = neZ d 

€k = h 2 k 2 /2m is the kinetic energy of one particle and a* k and a k are the usual boson creation and 
annihilation operators corresponding to momentum hk, satisfying the commutation relations 

[a k , a* k ,] = 5kk'i [ a k, o>kf\ — [ a L a *k'\ = 0- 
In both of our models the PBG (JIJ interacts with a one-mode photon field with Hamiltonian 

Qb*b, (2) 

where b* and b are the photon creation and annihilation operators, satisfying [b,b*] = 1, and 
Q > is the energy of a single photon. 

In both our models we assume that the photons interact with bosons linearly. The interac- 
tion in the first model has the form: 



whereas in the second one it is: 
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Hence the Hamiltonians of our models have the form : 



UiA = \gi{a* Q b + a Q b*), (3) 



U 2AL = ~g 2 {a* Q ,b* + a Q .b). (4) 



H 1:2A = f A + nb*b + U 1 , 2A , (5) 
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where T A := J^keA, k^o e k a *k a k and #1,2 > 0. 

The different forms of interactions (JHJ) and (J1J) come from two possible mechanisms of the 
light-bosons couplings. 

It is known that in the cubic box the conventional BEC of the PBG occurs only in the 
mode k = 0, see e.g. [10J. Since the scattering of the Q-mode light is most important on the 
macroscopic amount of condensed particles, following ^T] and ^2] we retain in the interaction 
only terms representing excitation from, and de-excitation back, to the BEC (in conformity 
with the minimal coupling in electrodynamics): 

Uia = (a* Q a b Q + a* a Q b* Q ), (6) 

with a coupling Ai > 0. Now assuming (as in [3]) that the depletion of the condensate can be 
neglected, one can simplify (jHJ) by substituting 

«o - vW> 4= a* - y/p^e^, (7) 



V w 

where po is the k = mode condensate density. This leads to the interaction (J3J) , with the cou- 
pling constant g\ = Xiy/po, after a gauge transformation eliminating the zero mode condensate 
phase if and after putting 6q = b. 

The interaction (j3J) has as its origin in the "four-wave mixing" mechanism , see e.g. pQ- 
P|. A condensate illuminated by a Q-mode "dressing" laser beam ("dressed condensate") can 
spontaneously emit pairs of photons and recoiling atoms. The simplest way to take this into 
account is described by the Hamiltonian 0: 

1 A 2 

U 2 \ = ( a Q' b Q" a obQ + a* b* Q a Q/ b Q „), (8) 

where 6q„, 6q» correspond to superradiated photons with a wave- vector Q" = Q — Q'. If as in jH] 
one neglects the depletion of the k = mode "dressed" condensate and the Q-mode "dressing" 
laser beam, then as in (j2J) the substitution of the corresponding operators by c-numbers gives 
instead of (jHJ) the interaction (jlj). Now g 2 is proportional to A 2v /po and the amplitude of the 
"dressing" laser beam, and 6q» = b. 

The aim of the present letter is to study the thermodynamic equilibrium properties of the 
models (0) and the possible phase transitions due to the coherent interaction of recoiled con- 
densate atoms with scattered (superradiated) photons. 

2. These models can be solved exactly by canonical transformations diagonalizing Hamiltoni- 
ans ©• To establish the thermodynamic properties we consider the grand- canonical ensemble 
Hamiltonians 

H lfiA (ji) = H li2A - t iN A , (9) 
where /x is the chemical potential and 

na= Yl G fc afc ' ( 10 ) 

fceA.fc^o 

is the particle number operator for the system excluding the ground state condensate. Since 
we shall follow a procedure analogous to that used for the PBG it is useful to recall some facts 
about the latter. 
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For the PBG the finite-volume grand-canonical thermodynamic functions exist only for 
/i < 0. Denoting the grand- canonical Gibbs state corresponding to the Hamiltonian Ta at 
inverse temperature (3 by (— ) T (/i), we have for p < 0, 



where N A = N A + a^a^ and 



1 f d d k 

pM = j^rL e^-i (12) 



For d > 3, poc = Po(0) < oo and therefore, to be able to access densities higher than p 0c in this 
case, one has to work with a given density rather than the chemical potential. We fix the finite 
volume chemical potential py by the equation 

" = (ir)<r>- (13) 

One then finds that for p < p 0c , p v f^oo where p^ < is the unique solution of the equation 

p = Po(n). (14) 

On the other hand for p > po c , 

^-mh^) +0(llv% (15) 

For the expectation of the zero mode occupation particle density we have 

Vm(p)w = \° f^;*' (16) 
y^oo \ V l Ta [p- poc for p > p 0c . 

2.1 Returning to our first model H 1A let us now define new boson operators and r]l,i]i 
by the canonical transformation: 

£i = cos-# aq + sini? b (17) 
771 = sin ■& oq — cos $ b . 

where 1? satisfies the equation 

tan 2# = - - — . (18) 

Q-e Q + p 

In terms of these operators we can write the Hamiltonian H 1A {p) in the form 

H 1A (p) = ^ k " ») a *k ak + + + ^iVi , (19) 



where 



E 1+ (p, n) = x - + t Q - p + ^(n + ^- eQ )2 + ^| 

Ei -{n, n) = l - jft + e Q - p - ^(Q + p - e Q ) 2 + g\ J . 



(20) 



Since Ei + (p, Q) > Ei-(p, Q), the thermodynamic stability of the Hamiltonian H\\(p) requires 
that 

E 1 .{p,Q)>0 or n<y. c (g x ,n) (21) 



where 

// c (^0)=min|o, -^ + e Q |. (22) 

Denoting the grand- canonical Gibbs state corresponding to the Hamiltonian H x \ by (—) HlA (//), 
we have for /i < /i c , 



V 2V [ \e^+ - 1 e^i- - I J E x +-E x - \e^+ - 1 e^- - 1 

(23) 

and 

for k 7^ 0, Q. Thus, since .&l_(/i, fl) > for < p c , we get for the total particle density 
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(/i) = ^^ ^fa-M)-i =^)- (25) 

Therefore as in the PBG we have to fix the finite volume chemical potential \iy by the equation 

P=(y) W. (26) 

\ / HiA 

and to distinguish two cases, p < p c and p > p c where 

pc = Po(jJ> c )- (27) 
Notice that by (|12|) if /i c < 0, then p c < oo even for d = 1, 2. 

(a) Case : p < p c . Using (|2l)j) we see that in this case p v —> Hoo, where p,^ < /i c is the unique 
solution of 

p = P o(p). (28) 

Since for these values of p, E x _(p, Q) is bounded below away from zero, by (J23j) we also get for 
the expectation of the Q mode occupation particle density 



Jim ( ) fa) = 0. (29) 

H 1A 



Similarly we get 



lim /^\ Qi V ) = (30) 

1 JY 1 | 1 V ^ + f 1 1 , , 

vToo IV \ \e^+ - 1 e^- - iy E 1+ - E x _ \e^+ - 1 eP*- - ' " 

for the density of photons. 
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(b) Case : p > p c . As above the analysis is similar to that for the PBG. Using (}2*3*|) - (|2T)j) we 

see that for V — - > oo, py takes the form: 

" A = "-F3(^) +0(1/n (31) 

and 

Ei-frv, n) = i 1 " + 0(1/V 2 ). (32) 

1/ ^(p - p c ) U-fJL c + t Q 

Then by (J23j) and (p?!^ we get instead of zero (see l2TJj) : 

Hm ( ) (p v ) = p - Pc7 (33) 

Hi A 



V^oo \ V 



implying the occurrence of BEC of the Bose gas in the Q-mode. Similarly, instead of zero as 
in (|5Uj) we obtain 

^(t)j^ = w^-^- (34) 

that is, one has condensation of photons occurring simultaneously with BEC (|33|). This corre- 
lation can also be seen from the limit of the boson-photon correlation (entangling) function: 

/a*b\ 

-0i(P _ Pc)/2ft for p>p c , ^ 
for p < p c . 



2.2 Now we condider our second model H 2 a We define new boson operators ££,£2 and 772,772 
by the canonical Bogoliubov transformation: 

£2 = cosh (ft aqi — sinh b (36) 
772 = cosh (j) b — sinh clqi . 



where satisfies the equation 



tanh 20 = - . (37) 

\l + 6q> - p 



In terms of these operators we can write the Hamiltonian H 2 a(p) in the form 
H 2 a(p) = ^2 ( e k~ n)a* k a k + E 2 + (p, fi)£a*6 + E 2 -(p, 0)775772 

+i{(£ 2 _+3 ! _)-( eg ,-, i + fi)} , (38) 
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where 



E 2+ {p, Q) = 1 | (e Q , - /i - fi) + ^( eQ ,-^ + fi)2_p2| (39 ) 

«) = 5 {-(eg' - A* ~ «) + \/ ( e Q' + g 2 } . 

The thermodynamic stability of the Hamiltonian H 2 \(p), always gives 

p<p c (g 2 ,n) (40) 



where 



p c (g 2 ,tt) =min{o, ~|^ + e Q'}- ( 41 ) 
p e = p(ji c ). (42) 



As in (J2ZJ) we let 

We have for p < p c , 

/ * V , x _ / / 1 , €Q/-// + fl / 1 1 

+ |i^-l M 43) 

£j 2 _|_ + -C/2 - 

(a*a*> Hj >) = e/3 ( efc - M) _ 1 for fc ^ 0, Q', (44) 
( & fe )ff 2 » = \ - ( el 3E 2+ _ I ~ e /3s 2 - _ i) + E 2+ + E 1 _ (e^+ _ i + ^s 2 _ _ i 

and 

<^>*p«) - 2^-,%-^ * - ir ' + ^ - ^ + o ■ (4e) 

If 6q/ > g 2 /4:Q then p c (g 2 ,Q) = and we return to the PBG with no condensation in the 
Q'-mode. If cq> < g 2 /AVL we have to study two cases g\ > AQ 2 and g\ < AQ 2 . Consider first 
the case g 2 > AVl 2 . 

(a) Case : p < p c Proceeding in a similar manner as for the first model we get 

lim /^9pL\ {flv) = i im /™\ M = Hm l°^P) (p v ) = 0. (47) 
y^oo \ V I H2A v^oo \ V I H2A v^oo \ V I H2A 

(b) Case : p > p c As for the first case here we obtain 

Km (/iy)=p _ Pc; (48) 
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V^oo \ V 



lim 




(49) 



and 



lim 




(50) 



In the case g\ < AVL 2 we find that the limits have the same form but with g 2 and Q interchanged. 
3. Clearly in our models the superradiance (see e.g. ,13 j and [H]) is directly related with 
BEC as is explicitly seen for example by comparing the formulae f)33|) and (|34p. Hence in spite 
of the simplicity of the models JHJ) they manifest an interesting cooperative phenomenon. The 
presence of the interaction between the Bose gas and radiation, compared to the PBG, occurs 
in both models at a lower critical density p c = po(eQ — g 2 /40) < poc = Po(0). Moreover the 
condensation in these models takes place not only in dimension d > 3, but also in dimensions 
d — 1 and d = 2. This shows clearly that the presence of radiation enhances the process of 
condensation in the Bose gas. 

It is also interesting to remark the value of the entangling boson-photon interaction energy 
in the two models, which one reads off from (J35)) and (J5U)) : 



For the first model based on the minimal coupling we obtain the negative interaction energy 
(bound state) in the presence of condensates, which is well-known [Zj-jH], [Hj- Whereas for 
the second model this interaction energy is positive, which is a completely different type of 
entanglement. 

These aspects of our results make contact with recent interests in entangled atom-photon 
states generated in BEC-superradiance experiments, see e.g. pQ, addressed to a variety of 
applications like tests of Bell inequalities, quantum cryptography and quantum teleportation. 
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